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It is shown that starting from one and the same transfer matrix formulation of QCD on a lattice, 
it is possible to obtain both the action of Hasenfratz and Karsch as well as an action where the 
chemical potential is not coupled to the temporal links. 



Numerical simulations of QCD at finite baryon density 
have met with long standing difficulties!! It is then nat- 
ural to inquire whether the situation cannot be improved 
by modifying the action used in the path integral. Special 
attention deserves the coupling of the chemical potential, 
whose form is due to Hasenfratz and Karschcl. They no- 
ticed that if in the euclidean path integral of relativistic 
gauge theories regularized on a lattice the chemical po- 
tential is coupled to the fermi fields as in classical sys- 
tems, specific counterterms are necessary which would 
be extremely unconvenient in numerical simulations. To 
avoid this shortcoming, guided by the analogy between 
chemical potential and temporal gauge held, they cou- 
pled the chemical potential in the same way as the gauge 
fields in the Wilson regularization. They showed that this 
definition does not require any counterterms, at least in 
the free theory. After some further analysis!, this defini- 
tion has been accepted as essentially unique for practical 
purposes!! 

The feature of the Hasenfratz-Karsh prescription 
which attracts our attention is the point splitting which 
implies that the chemical potential must be coupled also 
to the temporal gauge links. It has a topological con- 
sequence, in the sense that only the temporal Polyakov 
loops give contributions to the path integral which de- 
pend on the chemical potential. This looks like a mere 
lattice artifact, because it does not have any known coun- 
terpart in the continuum. But if the interplay between 
chemical potential and Polyakov loop is an artifact it 
could be eliminated. This appears an interesting pos- 
sibility from the point of view of numerical simulations 
because a term of toplogical nature is expected to affect 
them significantly. 

Therefore we started a systematic investigation of 
QCD at finite baryon density. We first avoided the use 
of the chemical potential by considering QCD in a given 
baryon sectord. The resulting effective action differs by 
the unconstrained action by terms with a quite differ- 
ent structure with respect to the coupling of the chem- 
ical potential in the Hasenfratz-Karsch action. Encour- 
aged by this result we reconsidered the definition of the 



chemical potential starting from the transfer matrix for- 
malism. The passage to the path integral, as far as the 
fcrmion fields are concerned, is based on the insertion 
of the identity expressed in terms of coherent states in 
the trace of the transfer matrix. But this can be done 
in different ways. A possibility which leads exactly to 
the Hasenfratz-Karsch prescription has been fully inves- 
tigatedB'El. The point splitting and the exponentiation 
of the chemical potential are in this case a consequence 
of the properties of Grassmannian kernels of fermionic 
operators!! But among other possibilities we show here 
that there is one where there is no point splitting and the 
chemical potential is not coupled to the temporal links. 
For Kogut-Susskind fermions in the flavour basis it is not 
coupled to the gauge fields at all. 

We start from the definiton of the grand canonical par- 
tition function at finite baryon density according to the 
ordered product 




ao is the lattice spacing in the temporal direction, uq 
the temporal component of the site position vector n, Nq 
is the number of links in the temporal direction, T the 
temperature 



Qb is the baryon charge operator and T q is the quark 
transfer matrix. The pure gauge part of the transfer ma- 
trix has been omitted because it does not play any role 
in the problem. The chemical potential \x is determined 
by the condition that the expectation value of the baryon 
number be 

q B = Z- l Tr | Q B expfaNo Qb) J[ %{n ) 1 . (3) 

T q can be writtenBi in terms of an auxiliary operator 
T q (n ) 
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T q (n Q ) = J-^^fqino + s), 



(4) 



where J is a function of the gauge fields which will be 
defined later and s = ±1 for Wilson/Kogut-Susskind 
fermions respectively. This difference in sign does not 
reflect any intrinsic difference, but is only due to the dif- 
ferent conventions adopted by Luschercl and Montvay- 
Miinsteiu, which we maintain for easy reference. The 
case of Kogut-Susskind fermions will be treated in the 
flavour basis, because we do not know a suitable for- 
mulation of the transfer matrix in the spin diagonal ba- 
sis. Obviously it is very simple to reexpress the result 
into this basis, but this introduces awkward nonminimal 
gauge couplings of the type discussed inEl. 

T q is defined in terms of quark-antiquark creation- 
annihilation operators x', y\ x, y acting in a Fock space. 
It depends on the time coordinate no only through the 
dependence on it of the gauge fields. In fact the creation 
and annihilation operators do not depend on no. They 
depend on the spatial coordinates n of the sites or, in the 
case of Kogut-Susskind in the flavour basis, of the blocks, 
and on Dirac, flavour and color indices, a, /, c sometimes 
comprehensively represented by I. 

In the transfer matrix formalism often one has to do 
with quantities at a given (euclidean) time no- For this 
reason we adopt a summation convention over spatial 
coordinates and intrinsic indices at fixed time. So for 
instance we will write 

i t M(n )x= ^ xl nI M m j. n j(n )x Ily j, (5) 

m,n,7, J 

where M is an arbitrary matrix. In this notation the 
baryonic charge Qb and the auxiliary operator T q (no) 
can be written 



y 1 y 



(0) 



T q (n ) = exp [-x t M(n ) x - y^M(n ) y] exp [y N(n Q ) x] . 

(7) 

Their form is the same for Wilson and Kogut-Susskind 
fermions but the matrices M and N are different in 
the two cases and will be specified later. The expres- 
sion of T 9 (no) is valid in the gauge Uq = A which is 
not admissible. But Menotti and Pelissettca extended 
the proof of reflection positivity to the gauge Uq ~ 11 
where Ua = 11 with the exception of a single time slice. 
In the construction of the path integral formulation of 
QCD at finite baryon density we do not need to fix the 
gauge, but to lighten the formalism we will nevertheless 
put Uq — 1 and we will reinstate Uq in the final result. 
The reader can check that keeping Uq in the intermedi- 
ate steps one arrives at the same result, provided some 
care is exercised: for instance when Uq ^ 1 the expres- 
sion x^M(no) x + y^M(no) y appearing in Eq. changes 
and does not commute with Qb any longer. We antici- 
pate that N is hermitean and also M is hermitean in the 
gauge U = 1. 



Before transforming the trace into a Berezin integral 
it is convenient to "distribute" the exponential of the 
charge operator among the factors in the trace 



Z = Tr l II ex P Q B 



(8) 



Now we introduce between the factors in the above equa- 
tion the identity 

l = J [dx + dx dy + dy] exp(— x + x — y + y)\xy >< xy\, 

(9) 



where the basis vectors 



\xy >= \ exp(— x x* — yy**) > 



(10) 



are coherent states and the x + ,x,y + ,y are Grassmann 
variables. They will be labeled by the time slice where 
the unit operator is introduced. For the other indices 
they are subject to the same convention as the creation 
and annihilation operators. 

Different realizations of the path integral arise depend- 
ing on the way this insertion is performed. The one lead- 
ing to the HaseHiratz-Karsch prescription has been dis- 
cussed in detailBQ. Here we consider another possibility 
based on the following way of writing the exponential of 
the charge operator 

exp(^ia Q) = J [dx + dx dy + dy] exp((55) 

• exp(— x + x — y + y)\x y >< xy\ (11) 

where 



SS = (1 — cosh(^oo)) (x + x + y + y) 
+ sinh(/iao)(x + x — y + y). 



(12) 



The above expression is obtained by expanding the ex- 
ponential of the charge operator, ordering all the anni- 
hilation operators to the left of the creation ones and 
inserting in each monomial the unity between the set of 
annihilation and the set of creation operators. For the 
rightmost exponential of the baryon charge one has to 
move the creation operators to the leftmost position be- 
fore taking the trace. Replacing the operators by their 
Grassmannian eigenvalues one gets the above result. 

Now the constojction of the path integral proceeds in 
the standard wayu, yielding the action for zero chemical 
potential plus the term generated by SS. To write down 
this term we must distinguish the case of Wilson fermions 
from that of Kogut-Susskind. 

For Wilson fermions we set s = 1 in Eq. [| and assume 



(13) 



M w (n ) =- In {[2KYB—1 (n ) 
N w (n ) = 2KB(n )-h(n Q )B(n Q ) 
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where K is the hopping parameter and 



B(no) = TL~KJ2 (u 3 M^ +) +T^U+(n Q )) 



3 = 1 



1 

c(no) = 5 X>i -lj^U+ (no)) ■ (14) 



3=1 



The link operators £7^ (no) have the standard Wilson 
variables U^(n) as spatial matrix elements 



V / m,n 



(15) 



We have introduced the translation operators Tj with 
matrix elements 



(2f>) 



3n2,ni±Cj 



(16) 



being the unit vectors 



(e,.)* = <W- (17) 
Next we define the quark field q by the transformation 



where 



x = B?P< ) +) q, y = B 1 -P ( Vl ] q 

q = ? f 7o 



P„ (±) = ^(l±7o). 



(18) 



(19) 



The jacobian of this transformation is the function J in- 
troduced in Eq. The partition function takes the form 

Zw = J [dqdq]exp(Sw + SSw) , (20) 

where 

S w = E {^E^fl + l„)U^n)q(n + e„) 

+g(n + e„)(l - 7/.)^(»)ff(»)] - ( 21 ) 

is the Wilson action with Wilson parameter r = 1 and 
zero chemical potential and the contribution of the chem- 
ical potential is obtained from Eq. |l2| 

5Sw = 2J {9(^0) [ (1 - cosh(//ao)) 

rto 

+ smh(fia )j ]B{n )q{n )}. (22) 

The first term breaks the chiral symmetry, but this is 
only a consequence of the breaking by the Wilson term. 
Indeed we will see that the corresponding term for Kogut- 
Susskind fermions respects the chiral invariance. Notice 



the "plus" sign in the exponential of the action, to com- 
ply with Liischer's convention. 

The case of Kogut-Susskind fermions will studied only 
in the flavour basis. The are now the block coordi- 
nates and the gauge fields are defined on the block links. 
We set s = — 1 in Eq. [| and assumeO 



M KS (n ) = 



N, 



<5 

Ks(n ) = { J2 [75 ® t 5 t 3 + 7j (p^Uiinoyr. 



3=1 



m 
K 



1® 1 + 75®*5<o}, (23) 



where m is the quark mass parameter and K the hopping 
parameter. In the tensor product the 7-matrices act on 
Dirac indices, while the ^-matrices 



(24) 



act on flavor indices. The projection operators P^ are 
given by 



P^ ±] = g I 1 ® 1 ± 7m7s (8 t B */*] • 
The quark field q is obtained by the transformation 
AVKP Q i+) q, y + = AVKP^q 



(25) 



x •■ 



(26) 



whose jacobian is the function J introduced in Eq. 
The partition function takes the form 



Z K s = J [dqdq] exp [-16 ( S KS + SS K s)} , 



(27) 



where 

Sks = t-^I ^2q( n ) (7m ® 1 - 75 ® *5*/i) U^{n)q(n + e M ) 

n 

+ 25(71)75 <8 ts*^ ?(") + amq(n)l ® ag(n)| (28) 

is the Kogut-Susskind action with zero chemical poten- 
tial and the chemical potential contribution is obtained 
from Eq. [TJ 

SSrs = --K"|y^g(n) [ (1 - cosh(/ido)) 75 ® £5*0 

n 

+ sinh(/ia )7o <S> l]g(n)|. (29) 

The factor 16 in front of the action accounts for the fact 
that the volume element with Kogut-Susskind fermions 
is 16 times larger than in the Wilson case. 

Let us emphasize that the partition functions with the 
present action and with that of Hasenfratz and Karsch 
are identical with one another. It is worth while checking 
this identity in the free case and comparing to the "naive" 
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definition. For simplicity we will consider only the Wil- 
son case neglecting the mass and the spatial Wilson term. 
But notice that we cannot omit also the temporal one, 
because otherwise we cannot construct the transfer ma- 
trix. We then consider the quark matrix 

Q = 7o(a V + a) + a al • V + -o§n + 1-R (30) 



where 



(V„/)(n) = — [/(n + e„)-/(r 



0V)(n) = — [/(n + e„) + /(n- 



0-2] 



(31) 



and 



(j = sinh(/iao), i? = cosh(/iao), new definition 

cr = /iao, i? = 1, "naive" definition. (32) 

We now evaluate the energy density £ at zero tempera- 
ture, setting et M = a. Normalizing at zero baryon density 



£ = - 



1 



2n 3 a 4 



d 6 qs 2 [1(h) -1(0)] 



(33) 



Therefore in the "naive" case the energy density has a 
quadratic divergence, while with the new definition 



(38) 



is equal to the value obtained by the Hasenfratz-Karsch 
prescription in the presence of the temporal Wilson term 
(without the spatial one). 

In conclusion we have shown how starting from one 
and the same transfer matrix we can realize identical 
path integrals with different actions for QCD at finite 
baryon density. In the action constructed in the present 
work the chemical potential is not coupled to the tem- 
poral links, in contrast to the Hasenfratz-Karsh prescrip- 
tion. It might therefore have quite different properties 
in numerical simulations which it could be interesting to 
investigate. 



where 



s 2 = C 1 - cos^j) . 



(34) 



1 

1(h) = — / dq [(smq - i<?) 2 + s 2 + (cos 9o - R) 2 ] ~ 



1 



- T [i-e(R-A)]. 



2\A4 2 - B 2 

6 is the step function. With the new definition 

R - A = ~ [s 2 + 2(1 - cosh(fia)] - i [s 2 - (/ia) 2 ] 



(35) 



A 2 - B 2 = 



1 , 
1 + tS z 



(36) 



while with the 'naive" definition 
1 



R-A=-[s 2 -(Ha) 2 ] 



A 2 -B 2 = s 2 



As 2 



(37) 
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